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CONVERSION  OF  EXIST1NC  CONTINUOUS  CONTROL  SYSTEMS  INTO  DIGITAL  CONTROL  SYSTEMS 


Kuldip  S.  Rattan 


Department  of  System  Engineering 
Wright  State  University 
Dayton,  OH  4543S 


ABSTRACT 


~f  The  problem  of  converting  existing  continuous- 
data  control  systems  into  digital  control  systems 
is  considered.  The  objective  of  this  paper  is  to 
develop  a  computer-aided  method  for  synthesizing 
the  pulse-transfer  function  of  the  digital  control¬ 
ler.  This  is  done  by  matching  the  frequency 
response  of  the  digical  control  system  Co  that  of 
the  continuous  system  with  a  minimum  weighted 
mean  square  error.  Formulas  for  computing  the 
parameters  of  the  digital  controller  are  obtained 
as  a  result.  The  design  technique  is  illustrated 
with  a  numerical  example  and  a  comparison  with 
previous  methods  is  also  presented.  <<£_. 

INTRODUCTION 

Digital  controls  are  the  subject  of  expanding 
Interest  in  both  the  academic  and  industrial 
communities.  This  is  especially  true  since  the 
advent  of  the  economical  microprocessor  and  its 
related  support  chips  which  have  made  digital 
controls  an  attractive  alternative  when  considering 
a  control  strategy.  For  the  past  few  years  new 
control  systems  have  been  designed  using  digital 
Instead  of  analog  controllers.  There  is  also  a 
need  for  converting  systems  which  were  designed  to 
be  controlled  by  continuous  controllers  into 
systems  to  be  controlled  by  digital  controllers. 
Conversion  of  existing  control  systems  is  a  problem 
of  considerable  Interest  to  various  branches  'of 
Industry. 

A  popular  method  for  converting  existing 
continuous  control  systems  into  digital  control 
systems  is  by  Tabak1.  This  method  converts  the 
transfer  function  of  the  continuous  controller 
into  the  z-transfer  function  of  the  digital  con¬ 
troller  using  prewar ped  bilinear  transformation. 

The  method  is  straightforward  and  easy  to  use. 

But  the  digital  control  system  obtained  by  this 
method  approximates  the  continuous  control  system 
only  when  the  sampling  frequency  is  sufficiently 
high  as  compared  with  the  highest  frequency  of  the 
continuous  system.  This  is  because  the  transfer 
function  of  the  plant  and  th^  feedback  are  not 
taken  into  consideration.  Thus  the  capabilities 
of  the  digital  controller  are  not  fully  utilized. 

Recently  the  author  of  this  paper  ant  Yeh2 
proposed  a  computer-aided  method  for  synthesizing 
•  digital  controller  by  matching  the  "discrete 


frequency  response"  of  the  digital  control  system 
to  the  frequency  response  of  the  continuous  model 
with  a  minimum  weighted  mean  square  error.  The 
"discrete  frequency  response”  of  the  digital  con¬ 
trol  system  is  obtained  by  putting  an  artificial 
sampler  of  duration  T  at  the  output  of  the  system 

and  then  substituting  z  •  e^“'r  In  the  overall 
pulse-transfer  function  of  the  system.  Even  though 
the  results  obtained  by  this  method  were  superior 
to  the  ones  obtained  by  Tabak's  method1,  this 
method  only  matched  the  response  of  the  systems  at 
the  sampling  instants.  The  purpose  of  this  study 
is  to  improve  on  this  method  by  matching  the 
"continuous  frequency  response"  of  the  digital 
control  system  with  that  of  the  continuous  model. 
The  "continuous  frequency  responses"  of  the  con¬ 
trol  systems  are  obtained  by  putting  an  artificial 
■ampler  of  duration  T>«0  at  the  output  and  Input 
of  the  systems,  substituting  in  the  overall 

transfer  functions  z  - - —  ,  s  -  w;  and  then 

1  “  T 

substituting  w  -  jy  (where  y  is  the  frequency  in 
the  w-doraain)  .  The  synethesis  of  the  digital 
controller  is  carried  out  by  minimizing  an  error 
between  the  "continuous  frequency  responses"  of 
the  digital  and  existing  systems  in  the  w-domaln. 

As  a  result  the  formulas  for  computing  the  para¬ 
meters  of  the  digital  controller  are  obtained. 

DIGITAL  CONTROLLER  DESIGN 
VIA  CONTINUOUS  FREQUENCY  MATCHING 

The  block  diagram  of  an  existing  continuous-  . 
data  control  system  may  be  drawn  as  shown  in 
Fig.  I. 


Fig.  1.  Existing  continuous-data  control 
system 

The  continuous  model  shown  in  Fig.  1  can  be  digl- 
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tallied  by  Inserting  a  sampler  of  duration  T  at  the 
error  input  e(t)  and  replacing  the  continuous  com¬ 
pensator  0-_  (s)  by  a  digital  controller  D(z)  and  a 

aero-order  hold  as  shown  In  Fig.  2.  The  sampling 

frequency  u ■  la  selected  to  be  sufficiently 

high  as  compared  to  the  highest  frequency  of  the 
continuous  system. 


response  of  the  continuous  model.  In  this  method 
we  only  matchad  the  response  et  the  sampling  In¬ 
stants.  It  Is  this  observation  which  prompted  the 
development  of  the  method  In  this  paper  which 
matches  the  "continuous  frequency  response"  of  both 
systems.  To  do  this  we  put  an  artificial  sampler 
of  duration  T^-0  secs  at  the  output  and  Input  of 

the  systems  shorn  in  Figures  1  and  2. 

Thus  the  overall  transfer  functions  of  the 
systems  shown  In  Figs.  1  and  2  with  a  sampler  T^-*0 

at  the  Input  and  output  can  be  written  as 
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Fig.  2.  Redesigned  digital  control  system 


The  overall  transfer  function  of  the  contin¬ 
uous  and  digital  systems  can  be  written  from  Figs. 
1  and  2  respectively  as 


Cm(s)  Cc(s)G(s) 

"rCO  "  F<s)  "  i+Gc<s)C(s)H(s) 

C.(s)  DT(s)G.  G(s) 

*4 -  -  C  (s)  - - - 

RT(s)  1+D(s)  G GHT(s) 


(1) 

(2) 


where  G  (s),  G.  (s),  G(s)  and  ll(s)  are  the  transfer 
c  n 

functions  of  the  continuous  controller,  zero-order 
hold,  plant  and  the  feedback  element  respectively. 
The  superscript  notation  denotes  the  sampling  rate 
T 

For  example,  R  indicates  that  the  signal  R  is 
sampled  at  ^  samples/second. 


It  Is  required  to  design  a  digital  controller 
so  that  the  steady  state  output  of  the  digital 
control  system  follows  che  desired  output  of  the 
continuous  model  for  all  sinusoidal  inputs  within 
the  frequency  range.  This  is  done  by  matching  the 
frequency  response  of  the  digital  control  system  to 
the  Ideal  frequency  response  of  the  continuous 
system. 


The  frequency  response  of  the  continuous  model 
can  be  obtained  by  substituting  s  ■  J«  in  eqn.  (1) 
and  varying  u  from  0  to  <*.  But  because  of  the 
hybrid  nature  of  eqn.  (2),  a  common  way  of  finding 
the  frequency  response  of  the  digital  control 
system  Is  to  put  an  artificial  sampler  of  sampling 
duration  T  secs,  at  the  output  of  the  system,  sub¬ 
stitute  z  "  in  the  overall  z-transfer  function 

and  then  vary  w  from  0  to  w/T.  We  will  call  this 
"discrete-frequency  response"  or  "sampled  spectrum". 
Recently  the  Author  and  Yeh2  proposed  a  method  for 
synthesizing  the  pulse-transfer  function  of  the 
digital  controller  by  matching  the  discrete  response 
of  the  digital  control  system  to  the  frequency 


Ti 

cd  (.) 

Tl 

R  (s) 


Ti 

C,  (.) 


T  T»  * 

«  (•)  C,C  (s) 


1-H>T(s)  GhCHT(s) 


<«> 


The  "continuous  frequency  response"  of  the 
continuous  and  digital  control  system  can  now  be 
obtained  by  transforming  eqns.  (3)  and  (4)  into  w 
domain  which  is  related  to  the  z -domain  by  a  bi¬ 
linear  transformation  (w  «  ^  and  the  s -domain 

by  w  -r  s  as  Tj-K)  (Whltback  and  Hofmann3)  and  then 
substituting  w  -  jy  in  the  overall  transfer  func¬ 
tions.  As  T x-*-0 ,  the  overall  transfer  functions  of 
the  continuous  and  digital  control  systems  given 
by  eqns.  (3)  and  (4)  can  be  written  In  the  v-domain 
as 


F(s)  -  F(s)  | 


(5) 


Cd(w) 


D(z)ChG(s) 


l+D(z)ChGH(z) 

S 


1+f 


(6) 


Equation  (6)  can  be  written  as 

C.(»)  --D(*>G(W>  - 

1  +  D(z)ChC!l(z) 


111 
1  "T 


(7) 


Let  the  transfer  function  of  the  digital  controller 
In  the  w-domaln  be  written  as 


D(w) 


av“+a  v™~l+ 

■I  m-i  i 

,  n . .  n-1. 

b  w  +b  ,w  + 
n  n-i 


♦•i***0 

1  » —  t  »<n 

4b^wfl 
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where  a's  and  b'a  arc  conacancs  to  be  determined. 
Substituting  v  •  j|  in  eqns.  (8)  and  (7),  we  get 


e(t)  -  F(Jt)-Cd(Jt) 


(15) 


D«»>  -  SSs  »> 

where 

2  16 

*  -  VV  +V  ~V  . . 

Q  -  •1-«jY2+*5Y*-*7Y*+ . 

t  -  -b2YJ'*^Y*-b6Y6+ . . 

M  -  b1-b3Y2+bJY4-b7Y6+  . 


F+1yQ 

*C1y)  - 


«l+JYi) 


1+iTL§m  <x*+iY*> 


(10) 


where 


C.CCJY)  -  C.C(8)|  -  X,+jY, 

a-)Y 

G^CHtjy)  -  GjOUz)  j  x  j  T/  "  X2+JY2 
a  -  1_JyT/2 


For  notatlonal  convenience,  let  the  numerator  and 
denominator  of  C^fjY)  be  denoted  by  H(y)  and  M(y) 
respectively,  l.e. 

Cd«Y>  '  SW  <U> 

In  terms  of  real  and  Imaginary  parts,  M(y)  and  N(y)’ 
can  be  written  as 


X(Y)  -  ♦  +  J*  (12) 

M(y)  -  o  +  Jt  (13) 

where 

♦  -  PX,-yOYi 

•  -  FYi+yVi 

o  -  1+L+PX2-Y<3Y2 
t  -  YM+PY2+YQX2 

The  transfer  function  F(Jy)  of  eqn.  (S)  can  also  be 
written  In  terms  of  real  and  imaginary  part  as 


An  optimal  set  of  a's  and  b's  can  therefore  be  ob¬ 
tained  by  minimizing  the  mean  square  of  the  magni¬ 
tude  of  c(y)>  Let  the  error  function  to  be  mini¬ 
mised  be  defined  as 

B  -  plFOY)  -  |2«*r  (16) 

,yl 


where  (yj'Y^)  is  the  bandwidth  in  the  v-domaln  on 

which  the  matching  is  required.  This  error  function 
is  to  be  minimized  by  adjusting  the  unknown  co¬ 
efficients  a.  and  b^.  However,  due  to. the  presence 
of  M(y)  in  the  denominator,  the  equations  obtained 
by  differentiating  E  with  respect  to  a1  and  b^  are 
nonlinear.  This  nonlinearity  may  be  removed 
by  modifying  the  error  function  (16)  as  . 

a* 

X  -  r2|F(jY>M(Y>  -  N(Y)|2dY  (17) 

This  integral  differs. from  that  in  eqn.  (16)  by  in- 

A 

eluding  |M(y) I  in  the  integrand  as  a  weighting 
factor.  This  modification  Is  not  new.  It  was 
suggested  by  Kalman4  as  an  identification  technique, 
was  used  by  Levy5  in  complex  curve  fitting,  Rao 
and  Lamba6  for  simplifying  linear  systems  and  by 
the  author  and  Yeh*  for  discretizing  continuous 
systems. 

The  integrand  in  eqn.  (17)  is  a  complex  func¬ 
tion  which  may  be  seperated  into  real  and  imaginary 

part  as 


|F(Jy)H(y)  -  N(y)|  -  Ky)  +  Jrn(Y) 

(18) 

where 

»(Y)  “  Ao-Bt-6 

(») 

■Cy)  ■  At+Bo-6 

(20) 

• 

By  substituting  eqns.  (18)-(20)  in  (17),  we 

*et 

E  -  [  2  [(Ao-Bt-6 )  2+  (At+Bo-6 )  2J  d  y 

(21) 

\ 

* 

Differentiating  with  respect  to  a's  and  b's  and 

equating  the  resulting  equations  equal  to  zero 

represents  (m+n+1)  linear  equations  In  (m+n+1) 

unknowns  a  to  a  and  b,  to  b_ .  These  equations 
o  ■  in 

'esn  be  written  in  the  matrix  fora  as  given  In  eqn. 

(22) 


F(JY)  -  A+Jl  (16) 

The  error  between  C.(Jy)  and  F(Jy) is  defined  ss 


Wh  -  J*2  (aV)  Yhdr 
1 

Eqns.  (22)-(26)  provide  a  design  algorithm  which 
has  been  programmed  on  a  digital  computer,  IBM 
(370/3031)  and  Is  available  in  the  form  of  a  sub¬ 
routine  written  In  FORTRAN  IV.  Note  that  the  order 
of  the  digital  controller  (integer  m  and  n)  is 
arbitrary.  Therefore,  using  this  algorithm,  one 
My  try  out  several  digital  controllers  of  dif¬ 
ferent  orders  and  select  the  best  design  by  con¬ 
sidering  the  tradeoff  between  cost  (complexity) 
and  performance. 


NUMERICAL  EXAMPLE 


The  example  considered  In  this  paper  is  to 
discretize  the  existing  continuous  control  system 
shown  In  Fig.  1.  The  transfer  function  of  the 
continuous  controller,  plant  and  the  feedback 
element  Is  given  by 


Cc(s) 


l40.J|16s 
HO.  139s 


The  magnitude  of  the  frequency  response  of  the 
continuous  model  is  obtained  using  eqn.  (30)  and  Is 
shown  in  Fig.  3.  Based  on  the  frequency  response, 
the  sampling  period  T  is  seteted  to  be  0.15s.  The 
transfer  function  of  the  zero-order  hold  is  given 
by 


Ch<"> 


s 


(27) 


(31) 
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Fig.  3.  Magnitude  of  the  frequency  response 
of  the  continuous  system 

From  eqns.  (28),  (29)  and  (31), 


GhGH(s) 


10(l-e~aT) 

»2(s+l) 


The  z- transform  associated  with  this  transfer  func¬ 
tion  is  given  by 


Substituting  w  ■  y  —•  in  eqn.  (37)  gives  the. 
pulse-transfer  function  of  the  digital  controller 


D(*)  • 


3.436t  -  2.191 
s  +  .2390 


The  frequency  responses  of  the  existing  contin¬ 
uous  system  and  the  digital  control  system  with  the 
digital  controller  given  by  eqn.  (38)  arc  shown  in 
Fig.  4-6.  It  can  be  seen  that  both  the  magnitude 
and  phase  angle  of  the  digital  system  matches 
closely  with  that  of  the  continuous  model.  If  the 
designer  is  not  satisfied  with  the  results  and 
wants  to  Improve  the  performance,  the  only  change 
which  needs  to  be  made  is  to  increase  the  order  of 
the  digital  controller.  The  algorithm  goes  through 
the  standard  steps  given  above  and  finds  the  para¬ 
meters  of  the  new  digital  controller.  Thus  the 
performance  may  be  improved  at  the  expense  of  in¬ 
creasing  complexity. 
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Substituting  r.  »  YTjyxT"  In  e<,n*  *33)*  we  8et 


Figure  4.  Nyquist  plot  of  the  frequency  res¬ 
ponses  of  the  continuous  system  and 
the  digital  control  system  designed 
by  different  methods. 


r  y  .  1Y  _  I.ISl-^w-.OOWw' 

chGH(jy)  -  x2+jy2  -  — (i+;99gj 


Substituting  s  ■  Jy  in  eqns  (28)  and  (30) ,  we  get 


C(JY>  -  X1+JY1  - 


w  -  Jy 


7«wOdwCwii»n  »rww  > 

►  •  •“  (Mf  Cora  oOr*  | 

*00  lefmd  Oar*  Co***om-  | *-1' 
l*vOr< Cnwrom* 

liMM  fivrimmi 
**■  O  CwMm 

Oeltr  timdum 


FOr)  -  A+JB  -  - 

wJ+8-194w*+37 • 142w+71-942 


•w  -  Jy 


Eqns.  (30),  (34)  -  (36)  are  the  only  information 
required  to  design  a  digital  controller  for  the 
digital  control  system  shown  in  Fig.  2.  Feeding 
these  eqns.  to  the  computer  program,  the  transfer 
function  of  the  first-order  digital  controller 
given  by  the  algorithm  in  the  w-doroain  is 


D(w)  - 


1.00516  +  . 34064w 


o  to  *o  to  so  100  120  no  iso  iso 

Frequency  in  rad/sec 

Figure  5.  Magnitude  plot  of  the  frequency 

responses  of  the  continuous  system 
and  the  digital  control  system 
designed  by  different  methods. 
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Figure  6.  Phase  plot  of  the  frequency  re¬ 
sponses  of  the  continuous  system 
and  the  digital  control  system  de¬ 
signed  by  different  methods. 


possible  with  that  of  the  continuous  model.  The 
psrsneters  of  the  digital  controller  are  obtained 
by  minimizing  a  modified  mean-square  of  the  error 
between  the  transfer  function  of  the  continuous 
and  the  digital  systems  in  the  w-domain.  An  ex¬ 
plicit  and  straightforward  futmula  tor  computing 
the  parameters  Is  obtained  and  has  been  programmed 
for  a  computer.  Since  the  design  procedure  Is 
computer-aided,  one  may  try  out  several  digital 
controllers  of  different  order  and  select*  the  best 
design  by  considering  the  trade-off  between  the 
coat  and  performance. 

The  numerical  example  considered  shows  that 
the  performance  of  the  digital  control  system  ob¬ 
tained  by  the  method  presented  is  superior  to  the 
ones  obtained  by  the  existing  methods. 
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COMPARISON  WITH  THE  EXISTING  METHODS 

The  digital  controller  for  the  digital  control 
system  discussed  in  the  numerical  example  is  re¬ 
designed  using  the  Tustin  transform7  (bilinear 
transformation) ;  methods  presented  by  Tabak1  (pre¬ 
warped  bilinear  transformation)  and  by  the  author 
of  the  paper  and  Yeh2  (discrete  frequency  matching). 
Using  these  methods,  the  continuous  controller 
given  In  eqn.  (27)  is  converted  into  digital  con¬ 
trollers  and  are  given  by 


Tustin  transform: 


D(z)  - 


2. 294z-1.5936 
*-.299 


Prewarped  bilinear  transform: 

,  2.42965z-l .68046 

D(I) - Z-0.-2509  (40) 

Using  the  method  presented  by  the  author  and 

Teh2: 

first  order: 

3.67872-2.31256 

»<*>  - 5083656  <41) 

'  second  order: 

.»(,)  ■  <«> 

x  -0.342-0.384 

The  frequency  responses  of  the  digital  control 
system  with  digital  controllers  given  by  eqns.  (39)- 
(42)  are  shown  In  Pigs.  4-6.  Comparison  of  the 
results  In  these  figures  demonstrate  that  the 
digital  controller  obtained  by  the  present  method 
gives  better  results. 

CONCLUSIONS 

A  method  for  converting  an  existing  continuous 
control  system  Into  a  digital  control  system  with 
similar  performances  Is  developed.  An  attempt  Is 
made  to  match  the  "continuous  frequency  response" 
of  the  digital  control  system  ns  closely  os 
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